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We present a comprehensive study of the energetics and magnetic properties of ZnO clusters doped 
with 3d transition metals (TM) using ab initio density functional calculations in the framework of 
generalized gradient approximation + Hubbard U (GGA+U) method. Our results within GGA+U 
for all 3d dopants except Ti indicate that antiferromagnetic interaction dominates in a neutral, 
defect-free cluster. Formation energies are calculated to identify the stable defects in the ZnO 
cluster. We have analyzed in details the role of these defects to stabilize ferromagnetism when the 
cluster is doped with Mn, Fe, and Co. Our calculations reveal that in the presence of charged defects 
the transition metal atoms residing at the surface of the cluster may have an unusual oxidation 
state, that plays an important role to render the cluster ferromagnetic. Defect induced magnetism 
in ZnO clusters without any TM dopants is also analyzed. These results on ZnO clusters may 
have significant contributions in the nanoengineering of defects to achieve desired ferromagnetic 
properties for spintronic applications. 
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I. INTRODUCTION 

Functional materials in reduced dimension are of im- 
mense interest from the point of view of both fundamen- 
tal science and technological applications as novel prop- 
erties emerge due to the complex quantum mechanical 
interactions at this size regime. One of the major appli- 
cations is in high density storage in miniaturized devices. 
Quantum confinement plays a key role in tailoring prop- 
erties of these small systems, e.g., atomic clusters con- 
taining tens of atoms. In this context, magnetic clusters 
and nanosized systems, either in a pure atomic assembly 
or as magnetically doped semiconducting systems, have 
attracted a lot of attention. An exciting material is tran- 
sition metal (TM) doped ZnO nanocrystals where fer- 
romagnetism is reported even at high temperatures;^^— 
This system has received special attention as naturally 
abundant environment friendly ZnO with its wide band 
gap coupled with magnetism may provide a unique op- 
portunity to achieve electronic, magneto-electronic, and 
opto-electronic multifunctionality in a single system. 

First principles electronic structure calculations on TM 
doped ZnO clusters have been performed to provide use- 
ful insights regarding the electronic structure and mag- 
netic interactions. Reber et al. have studied a small 
cluster of ZnO (Zni20i2) in the framework of general- 
ized gradient approximation (GGA) and reported that 
doping of Co leads to ferromagnetic coupling due to di- 
rect exchange interaction in this system^. It was also re- 
ported that doping of Mn into the Zni20i2 cluster yields 
antiferromagnetic (AFM) ground state for small Mn-Mn 
distance while ferromagnetic and AFM states are degen- 
erate for large Mn-Mn distance^ . Very recently the struc- 
tural and magnetic properties of several 3d TM doped 
Zni20i2 clusters were investigated in the framework of 



GGA.^ The authors find that most of the TM dopants 
stabilize either in paramagnetic or in antiferromagnetic 
state, except Co and Ni, where the ferromagnetic and an- 
tiferromagnetic states are found to be energetically very 
close. 

The experimental works on TM doped ZnO clusters 
are primarily focused on the study of ferromagnetism 
(FM) in TM doped ZnO nanoparticles.^"*' Recently Kar- 
makar et al. have found ferromagnetism in Fe doped 
ZnO nano-crystalsi^ However for these systems the local 
probes like electron paramagnetic resonance (EPR) and 
electron spin resonance (ESR) reveal an unusual ionic 
state of the dopant Fe atoms. They find the presence of 
both: Fe^+ and Fe"^+ ionic states in the system and argue 
that these are crucial for ferromagnetism. The presence 
of Zn vacancies are speculated to play an important role 
in stabilizing Fe'^+ ions in the system. X-ray Magnetic 
Circular Dichroism (XMCD) studies on the same samples 
later confirmed that Fe'^"'" ions are located at the surface 
of the nanoparticles and are crucial for their magnetic 
properties^ It has been reported that the type of car- 
riers also play an important role in ferromagnetism in 
TM doped ZnO quantum dots: Mn^"*" doped ZnO quan- 
tum dots become ferromagnetic in presence of p-type 
carriers, while n-type carriers induce ferromagnetism in 
Co^^ doped ZnO quantum dots.'^ Interestingly an ex- 
perimental report indicates that ferromagnetism in ZnO 
and in many other oxide nanoparticles may be realized 
even in the absence of TM doping presumably promoted 
by defectsiiSi The possibility of inducing room temper- 
ature ferromagnetic-like behavior in ZnO nanoparticles 
without doping magnetic impurities but coating with N 
and S containing ligands such as dodecylamine and do- 
decanethiol has also been realizedi^i 

The above discussion clearly suggests that magnetism 
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in TM doped ZnO clusters and nanosystems is rather 
subtle, where surface effects and defects play a crucial 
role to stabilize FM. Therefore, an understanding of the 
defect induced magnetism is crucial to design these ma- 
terials for desirable property. In a recent work using 
GGA+U on Fe doped ZnO clusters we have shown that 
defects under suitable conditions can induce ferromag- 
netic interactions between the dopant Fe atoms, whereas 
antiferromagnetic coupling dominates in a neutral defect- 
free cluster In addition, there are several calculations 
that suggest that ZnO nanocluster may be ferromagnetic 
without TM dopants. Schoenhalz et al. have proposed 
that the surface states promoted by extended defects {e.g. 
grain boundaries) may play an important role in mediat- 
ing ferromagnetic interaction in such materials.— A re- 
cent first principles study of ZnO nanostructures coated 
with ligands containing N and S also reported ferromag- 
netism due to the redistribution of charge promoted by 
ligand capping. 

In this paper we shall present a systematic study and 
analysis of the electronic structure and magnetism of 
transition metal doped Zni20i2 cluster in the framework 
of GGA-I-U. We have carried out our studies both in the 
presence as well as in the absence of defects in order 
to clarify the nature and origin of magnetic interactions. 
Our study along the id transition metal series will be use- 
ful to extract trends. Our calculations demonstrate that 
doping TM except Ti in pristine ZnO does not render 
the cluster ferromagnetic, however defects under suitable 
condition can induce ferromagnetic interaction between 
the dopant TM ions. 

The remaining part of the paper is organized as follows: 
In section |lT] we have described the technical details and 
the method of our calculations. Details about obtain- 
ing the structure of the cluster is given in section Hill We 
have discussed the results for doping of TM into the clus- 
ter in section lTV Al Effect of defects namely Zn and O va- 
cancies with different charge states in the presence and in 
the absence of TM dopants are discussed in sections HVBI 
and |IVC] respectively. Finally, we have summarized our 
work in section |Vl 



II. METHODOLOGY 

The total energy and electronic structure calculations 
presented in this paper are carried out in the frame- 
work of ab initio density functional theory (DFT). We 
have used plane wave basis along with the projector aug- 
mented wave (PAW) method^^ as implemented in the 
VASP— code. The exchange-correlation term in DFT 
was treated within the generalized gradient approxima- 
tion (GGA) due to Perdew-Burke-Ernzerhof (PBE).!^ 
For the description of dopant TM atoms we have em- 
ployed GGA-I-U with the typical value of on-site d — d 
Coulomb interaction U = 4.0 eV and the onsite ex- 
change interaction J = 1.0 eV. PAW potentials with 
12 valence electrons {3d^^4s^) for Zn, 6 {2p^2s^) for O, 



and 2 + n (3d"4s2); n = 2 to 9 for TMs (from Ti to 
Cu) have been used. The plane wave energy cutoff was 
taken to be 500 eV. The cluster was simulated in a large 
cell (a cube of volume (20 A)'^) to ensure sufficient sep- 
aration between the periodic images. Only one /c-point 
{i.e. r point) was used for these calculations. We have 
relaxed the atomic positions to minimize the Hellman- 
Feynman forces on each atom with the tolerance value of 
0.01 eV/A. 

The formation energy (FE) of defects is the energy 
for the reaction needed to create the defect from the 
ideal material. It is calculated using the following 
expression"'^^: 

FE - £;[Zn„,0„(a,g)] -S[ZnpOp(pure)] 

+ ^UafJ-a + q{Ey + ep), (1) 
a 

where a is the defect atom added to or removed from 
the pristine ZnO cluster, is the number of each de- 
fect atom: = —1 (-f-1) for adding (removing) one 
atom, E[Zii^On{a,q)] is the total energy with defect a 
and charge q, while £'[ZnpOp(pure)] is the total energy 
of pure ZnO cluster, e^? is the Fermi level measured with 
respect to the energy {Ey) of the highest occupied molec- 
ular orbital (HOMO) in pristine ZnO cluster, fia is the 
chemical potential of atom a in some suitable reservoir. 

The chemical potential of a particular atom is the en- 
ergy supplied (released) to add (remove) one atom of that 
type to (from) the system. It depends on the conditions 
under which the material has been prepared. However, 
there are limits on the values of chemical potentials of Zn 
and O for the formation of ZnO cluster. In equilibrium 
the chemical potentials of all Zn and O atoms should add 
up to the total energy of the cluster, i.e. the following 
equation must be obeyed: 

p{^J■Zn + fJ-o) = £^[ZnpOp(pure)], (2) 

where p is the number of ZnO units in the cluster. How- 
ever, the chemical potential of an individual atom can 
not be more than its standard value in elemental struc- 
ture because otherwise the standard elemental structure 
will be more stable and in that case the cluster will never 
be formed. This imposes the following constraints on the 
chemical potentials: 

Mzn < Mzn and Mo < Mo' (3) 

where and /Iq are chemical potentials in standard 
elemental structures of Zn and O respectively. Both the 
above constraints can be simultaneously satisfied if fi^^ + 
Mo — '^[ZnpOp(pure)]/p, and Eq. ^ can be rewritten in 
the following way: 

P(/4n + Aiin + Aio + /^o) = -B[ZnpOp(pure)], (4) 

where is the excess chemical potential of Zn, defined 
as: Mzn = ^^Zn — t^'zn- Excess chemical potential of O 
(Mo) is also defined in the similar way. We gather from 





Eq. ([3]) that the conditions: < and iiq < must 
be satisfied and Eq. Q indicates that there is only one 
free parameter between fi^^ and /Kq there are lower 
limits on the values of these excess chemical potentials. 
Here we have taken two extreme limits of chemical poten- 
tials to calculate the formation energies: (i) Zn rich limit 
= 0, fiQ is minimum) and (ii) O rich limit {fi^ = 0, 
is minimum). The chemical potential in standard 
elemental structure of atom a {fi% ) was taken to be the 
energy per atom of the element in the most stable ele- 
mental structure (bulk solid or molecular form) . We have 
simulated different defect charge states by adding or re- 
moving electrons from the system. The additional charge 
was compensated by a uniform jellium background suit- 
able for the system to ensure charge neutrality of the 
system. In view of the large supercell chosen in our cal- 
culation, we have not included any finite size correction 
for the charged systems 




III. STRUCTURE OF THE CLUSTER 



FIG. 1: Optimized structure of Zni20i2 cluster. 



We have considered Zni20i2 cluster for our simula- 
tions. The cluster of this particular size has been re- 
ported to be the magic size cluster with the most sta- 
ble configuration<ii2£ In this context, it is interesting to 
note that the experimental absorption energies of small 
ZnO nanoparticles are found to be very close to the op- 
tical excitation energy of Zni20i2 cluster calculated in 
the framework of time dependent density functional the- 
ory (TpDFT)Jl This observation leads the authors of 
Ref. [2l| to conjecture that the surface structure of small 
ZnO nanoparticles is very similar to that of Zni2 0i2 clus- 
ter. In view of the above discussion, our study of TM 
doped Zni20i2 cluster may help understanding the role 
of surface in the magnetic properties of TM doped ZnO 
nanoparticles. 

We have considered 12 ZnO units in several different 
initial configurations and relaxed the atomic positions 
to obtain the ground state structure by minimizing the 
forces on the atoms. The lowest energy structure con- 
taining all the atoms arranged on the surface of a hollow 
sphere is displayed in Fig. 1. Each atom of the cluster 
is coordinated with three neighbors and the atoms are 
arranged periodically with 6 four-atom rings and 8 six- 
atom rings. This structure is in good agreement with 
the structure reported earlier for the same size of the 
cluster— The diameter of the cluster is ^^6.35 A. 
Although the size of the cluster is much smaller than the 
nanoparticles synthesized experimentally, however since 
all the atoms are located on the surface of a sphere in 
this cluster, it provides a unique opportunity to study 
in details the surface contribution to magnetism that 
will be relevant for a ZnO nanoparticle. The gap be- 
tween HOMO and the lowest unoccupied molecular or- 
bital (LUMO) is calculated to be -2.34 eV. 



IV. RESULTS AND DISCUSSIONS 

As a next step, we consider doping of the cluster with 
TM ions. In this context, recent theoretical studies on 
TM doped bulk ZnO indicate that standard exchange- 
correlation (XC) functionals such as local spin density 
approximation (LSDA) /GGA are in general not adequate 
to describe these systems. One such report^ points out 
that the LSDA predictions might be misleading as the 
localized character of TM impurities are not taken into 
account. On the contrary, incorporation of Hubbard U 
gives more appropriate description of the system~. Very 
recently, Iu§an et al<^ have shown that the magnetic 
exchange parameters calculated with LSDA+U for Co 
doped ZnO are in good agreement with the magnetic 
measurementsi^ Considering the importance of Coulomb 
correlation, we have employed GGA+U method for the 
description of TM-d states in the TM doped ZnO cluster. 
We have calculated the stability of the cluster on doping, 
the corresponding magnetic moments, and magnetic in- 
teractions between two TM ions. The positions of the 
atoms in the cluster were fully relaxed upon doping. 



A. Doping of TM: GGA+U results 

We have only simulated substitutional doping in view 
of a recent report on Mn doped ZnO cluster where it is 
found that substitutional doping of Mn at the Zn site is 
energetically more favorable compared to endohedral or 
exohedral dopingSi. From our calculations we find that 
doping of TM at the Zn site is energetically more favor- 
able than doping it at the O site. Therefore we have 
considered only substitutional doping of TM at Zn site 
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Ti V Cr Mn Fe Co Ni Cu 



FIG. 2: Formation energies (in Zn rich and O rich limit) and 
magnetic moments are shown for Zni2 0i2 cluster doped with 
diflerent TMs. 

for our work. 

We have calculated the formation energies of doping 
one TM ion into the system using Eq. (IT|) both at Zn 
rich and O rich limits. Magnetic moments correspond- 
ing to doping of different TMs are also calculated. The 
results of our calculations are presented in Fig. 2. We 
find that the Zn rich limit is hardly favorable for TM 
doping. However, as expected, at O rich limit formation 
energies corresponding to doping of all TM ions are neg- 
ative, which implies that O rich limit is more favorable 
for doping of TM into the cluster. In practice, doping of 
Ti in ZnO at the O rich limit may not be achieved easily 
as there is a possibility of formation of other oxides of Ti 
in this limit. 

The magnetic moments of the clusters show interest- 
ing features depending upon the dopant TM ions. As 
we move from Ti to Cu, the magnetic moment increases 
first, becomes maximum (5 /is) for Mn and then de- 
creases monotonically upto Cu. The magnetic moments 
are mainly localized on the TM ions. The TM ions are 
expected to be in 2-1- oxidation state since they substi- 
tute Zn^+ ion. The values of magnetic moments can be 
explained by Hund's rule i.e. assuming high spin arrange- 
ment of 3d electrons for TM^+ ions, with the exception of 
Ti. The magnetic moment is completely quenched upon 
Ti doping. It is interesting to note that, as opposed to 
other TM dopants, the Ti ion and the three neighboring 
oxygens surrounding Ti are confined in a plane upon re- 
laxation. The Ti-0 bond lengths are found to be quite 
small (1.86, 1.86, and 1.89 A) compared to other TM-0 
bond lengths {e.g. Mn-0 bond lengths are 1.95, 2.04, 
and 2.04 A). Therefore, the crystal field experienced by 
the Ti ion is stronger compared to that for the other TM 
dopants. This crystal field is in fact stronger than the 
exchange splitting for Ti, thus stabilizing the low spin 
state for the singly doped Ti ion. 

The spin polarized densities of states corresponding to 
doping of one TM ion into the ZnO cluster are shown 
in Fig. 3. The DOS corresponding to the pristine ZnO 
cluster has also been shown at the rightmost panel for 
comparison. The figure shows that the TM-d states are 
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FIG. 3: Total density of states (DOS) and partial density 
of states (PDOS) corresponding to TM-d states for Zni20i2 
clusters doped with one TM. Partial Density of States of TM- 
d states are shaded, whereas solid lines indicate total density 
of states for ZnO cluster doped with one TM. For each verti- 
cal panel, sub-panel of left hand side indicates majority spin 
states and sub-panel of right hand side indicates minority spin 
states. DOS of pristine ZnO cluster has also been shown at the 
rightmost panel for comparison. Gaussian smearing scheme 
with 0.1 eV of smearing width was used for the DOS calcula- 
tions. 



lying deep into the gap region. We gather from the figure 
that as expected for high spin TM^+ ions, the minority 
(i-states are completely unoccupied for V, Cr, and Mn, 
while for Fe, Co, Ni, and Cu the minority states are pro- 
gressively filled, accounting for the reduction in magnetic 
moment. Further, the overlap of the TM-d states with 
the host states increase as we move from Ti to Cu. 

After exploring the possibility of doping of TM into 
the cluster, we have calculated the magnetic exchange 
interaction between TM ions in this system. To study 
the exchange interaction between TM ions it is essential 
to substitute at least two TM ions into the cluster. We 
have identified 7 different possible distances between Zn 
atoms, where a pair of TMs can be doped in pristine 
ZnO. Hence we have got 7 different configurations with 
increasing distances. Doping of TM ions in all possi- 
ble configurations are assumed and the atomic positions 
are relaxed separately for parallel and antiparallel orien- 
tations of magnetic moments on TM ions. The energy 
differences between antiparallel (AFM) and parallel (fer- 
romagnetic) states for different configurations have been 
plotted as a function of transition metal and displayed 
in Fig. 4. Since for higher configuration numbers (larger 
distances between TM ions) the exchange interactions 
effectively die out, here we have only shown the interac- 
tions upto fourth configuration. The fast decay of the 
interaction strength with distance may be attributed to 
the following reason: It has been observed that TM-d 
states are mainly located deep inside the HOMO-LUMO 
gap region (see Fig. 3). Since these states are not de- 
localized over the cluster, they can only interact with a 
state located very close to it in space, accounting for the 
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FIG. 4: Exchange interaction as a function of dopant TM for 
first four configurations. 



FIG. 5: Partial DOS corresponding to Ti-d states near the 
Fermi level for Ti2ZnioOi2 cluster. 



sharp decay. 

We find that except Ti, the dominant exchange inter- 
action for all the other transition elements (V-Cu) doped 
into ZnO cluster to be antiferromagnetic. This is in sharp 
contrast to the GGA results discussed in the literature ji"— 
where one finds the tendency of ferromagnetism in Co, 
Cu, and Ni doped ZnO cluster and the system is non- 
magnetic upon Ti doping. Therefore our calculation em- 
phasizes the importance of inclusion of Hubbard U in the 
study of magnetism in TM doped ZnO cluster. 

Fig. 4 reveals that the exchange interaction for a pair 
of Ti atoms is ferromagnetic. It has a magnetic moment 
of 2 fj,B- We have seen that Ti has a large crystal field 
splitting and weak exchange splitting. However, a pair of 
Ti atoms may interact forming bonding and antibonding 
states and this splitting is rather strong for Ti in compar- 
ison to other TM atoms doped into the system. In con- 
trast to other TM's, a pair of Ti atoms are found to come 
closer to each other upon relaxation. For example, the 
separation between Ti atoms for the first three configura- 
tions are 2.44, 2.66, and 2.75 Arespectively, while for Mn 
it is 2.91, 3.37, and 3.38 A. This provides strong bonding- 
antibonding interaction between similar spin states pro- 
moting an arrangement of majority and minority states 
as shown in Fig. 5 (partial DOS corresponding to Ti-c? 
states), accounting for the calculated magnetic moment. 



B. Effect of vacancy 

From the results of the preceding section we gather 
that ferromagnetism is very unlikely to be stable in tran- 
sition metal doped pristine ZnO clusters with the possible 
exception of Ti. However, there are several reports that 
claim room temperature ferromagnetism in Mn, Fe and 
Co doped ZnO, which suggests that defects may play an 
important role in stabilizing ferromagnetism in such sys- 
tems. Different types of native defects {e.g. interstitials, 
vacancies, and antisites) have been identified for ZnO in 
bulk form. It was reported from first-principles calcu- 



lations that Zn and O vacancies are the most relevant 
defects in ZnO^S. 

In view of the above we have studied the TM doped 
Zni20i2 clusters in the presence of defects. Here we have 
only considered Mn, Fe, and Co as dopants, since major- 
ity of the experimental results on ZnO nanosystems are 
reported with these dopants. 

To find out the stable defect states in TM doped ZnO 
cluster, we study the formation energies for a single tran- 
sition metal doping with and without Zn and O vacancy 
in different charge states. Similar to the previous cases, 
we have calculated formation energies in Zn and O rich 
limit. We have studied two different configurations de- 
pending on the distance between TM ion and vacancy, 
namely near (TM ion and vacancy are at the closest pos- 
sible positions) and far (distance between TM ion and 
vacancy is maximum). We considered both Zn rich and 
O rich conditions. The Fermi level has been varied from 
the HOMO level of the pristine ZnO cluster to the value 
of the experimentally reported band gap of bulk ZnO. 

The formation energies for doping of Mn, Fe, and Co 
into the ZnO cluster with and without vacancies at differ- 
ent charge states are shown in Fig. 6. Since the formation 
energies in the Zn rich limit is never found to be lower 
than that in the O rich limit, here we have only displayed 
our results in the latter limit. Fig. 6 shows that TM 
doped ZnO in the charge state -1-1 has the lowest forma- 
tion energy at the p-type regime, whereas Zn vacancy in 
the charge state —1 is the most stable state at the n-type 
regime. Here we note that the stable defects may lead to 
3+ oxidation state of the TM atom on the surface of the 
cluster. This is in agreement with a recent experimen- 
tal report on Fe doped ZnO nanocrystals^. The unusual 
ionic state may be attributed to different coordination of 
TM atoms on the surface of the cluster compared to the 
bulk terminated ones due to the loss of ligands. 

Having established the most stable defects in the TM 
doped cluster, we have investigated whether the most 
stable configuration may support the tendency of FM. 
The energy differences between the AFM and ferromag- 
netic configurations with the inclusion of two TM ions in 
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FIG. 6: (Color online) Calculated formation energies for sin- 
gle TM atom doped cluster at oxygen rich limit with and 
without vacancies in different charge states. Only the stable 
part of a particular charged state is shown. When vacancy 
is present, TM atom and vacancy are at the closest possi- 
ble positions. The Fermi energy ep has been varied upto 
experimental band gap of bulk ZnO and dotted vertical line 
indicates the calculated band gap for pure cluster. 



the absence and presence of Zn and O vacancies in vari- 
ous charge states are shown in Fig. 7. Here the TM ions 
are always assumed to be closest to each other. Calcula- 
tions are carried out when the defects (Zn or O vacancy) 
are near as well as far from the pair of transition metals 
(TM) doped into the cluster. From Fig. 7 we gather that 
in the absence of defects Mn, Fe, Co doped ZnO cluster in 
the charge state exhibit tendency for ferromagnetism. 
Zn vacancies in the —1 charge state also induces a ten- 
dency of FM into the cluster for all three dopant transi- 
tion elements (Mn, Fe, and Co). Although Zn vacancy in 
charge state shows tendency of ferromagnetism upon 
doping of Mn, this result is not interesting as Zn vacancy 
in -1-1 charge state is not stable. We note that in either 
case, TM doped ZnO cluster in the charge state with- 
out any defects or in the charge state -1 in the presence 
of Zn vacancy leads to 3-|- oxidation state for the dopant 
TM and are important to promote ferromagnetism. 

In order to understand the possible mechanism that in- 
duces ferromagnetic order, we consider a representative 
case namely a cluster doped with 2 Mn ions along with 
one Zn vacancy in the —1 charge state. The presence of 
a Zn vacancy in the Zni20i2 cluster doped with a pair of 
Mn atoms in the charge state -1 relaxes significantly, as 
shown in Fig. 8(a). In the present case the oxygen atoms 
surrounding the Zn vacancy are found to undergo large 
(^^28%) outward relaxation. Similar trend is seen for 
the other stable charge states. Next we shall address the 
magnetism of the cluster. In the presence of a Zn vacancy 
in the charge state —1 there is one hole in the system. 
When Mn ions are substitutionally doped at Zn sites, its 
oxidation state should be Mn^+ . Due to deficiency of one 
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FIG. 7: Calculated energy differences between AFM and FM 
orientations for different charge states. Positive (negative) 
value indicates ferromagnetic (antiferromagnetic) ground 
state. 
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FIG. 8: ZnO cluster doped with 2 Mn ions along with one 
Zn vacancy and —1 charged state, (a): relaxed structure of 
ZnO cluster doped with two Mn atoms and one zinc vacancy 
in the charge state -1, (b): level diagram showing expected 
interaction between d states of two Mn atoms (energy levels in 
black correspond to majority spin states and energy levels in 
red correspond to minority spin states), and (c): partial DOS 
for Mn-d states (inset shows the charge density corresponding 
to highest occupied Mn-d (majority) state). 
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electron in the system, one Mn ion may be in 3+ oxida- 
tion state resulting in mixed valency with a combination 
of Mn^"*" and Mn'^+ states present in the system. In such 
a situation if the spins are ferromagnetically aligned as 
shown in Fig. 8(b) then there is a possibility of lower- 
ing energy by spin conserved hopping as illustrated in 
Fig. 8(b). To verify the above mentioned mechanism, we 
have plotted partial DOS for Mn-d states in Fig. 8(c) and 
the charge density corresponding to the highest occupied 
Mn-d (majority) state in Fig. 8(c)(inset). We gather from 
Fig. 8(c) that the minority Mn 3d states are completely 
unoccupied, while the majority states are completely oc- 
cupied except for one state, resulting in a net moment of 
9 /LtB and consistent with the mixed valent states as shown 
schematically in Fig. 8(b). A plot of the charge density 
in a small energy window corresponding to a single Mn-d 
state in the majority spin channel (Fig. 8(c)(inset)) re- 
veal that both the Mn ions, as well as the nearby O ions 
contribute to this state, indicating that the electron is 
hopping between both the Mn ions either directly or via 
O ions. Thus hopping induced interaction between the 
dopant Mn atoms stabilizes the tendency of FM for this 
system, similar to our previous observation for Fe doped 
ZnO cluster with the same defect state^'^. 

In this context, an experimental report has also 
pointed out that deficiency of electrons can induce fer- 
romagnetism in Mn doped ZnO quantum dots."^ A calcu- 
lation by Feng et. al. in the framework of TDDFT also 
claims that double exchange mechanism induced by opti- 
cal excitation can stabilize ferromagnetism in Mn doped 
ZnO quantum dotsi^ In their description, optical exci- 
tation generates holes in the valence band of the Mn ion, 
thus making it behave like Mn^+ ion. Holes thus gen- 
erated participate in the double exchange process and 
thereby stabilize ferromagnetism. The above discussion 
points to the fact that the coexistence of Mn^+ and Mn'^+ 
ions in the system may be important for ferromagnetism. 
In the present work, we have demonstrated the role of de- 
fects to satisfy the requirement for mixed valency of Mn 
in the system. 

Although oxygen vacancies are not very stable in the 
system, it has interesting effects on magnetic interac- 
tions. This is particularly visible (Fig. 7) in case of Co 
doping where both -1-1 and —1 charges induce ferromag- 
netic coupling between Co atoms. Once again, these re- 
sults are similar to a recent theoretical calculation on the 
feromagnetism induced by O vacancies in a Co doped 
ZnO system3^. 



C. Defects without TM doping 

So far we have studied the role of defects in influ- 
encing the magnetic interactions between the dopant 
TM ions. However, recent experimental reports on ZnO 
nanoparticlesiS and ZnO thin films^^ that ZnO may be 
ferromagnetic even in the absence of TM dopants moti- 
vated us to investigate if defects {viz. Zn and O vacancies 
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FIG. 9: (Color online) Calculated formation energies for Zn 
and O vacancies in the pristine ZnO cluster at Zn rich and 
O rich limit. Only the stable part of a particular charged 
state is shown. The Fermi energy ep has been varied upto 
the experimental band gap of bulk ZnO and the dotted ver- 
tical line indicates the calculated HOMO-LUMO gap of the 
pristine ZnO cluster. 



in different charge states) can promote ferromagnetism in 
ZnO cluster. For this purpose, we have calculated the for- 
mation energies corresponding to these defects and the 
results are displayed in Fig. 9. We gather from this fig- 
ure that O vacancies in the Zn-rich limit are stable in 
the charged state and +1, while the Zn vacancies in 
the O rich limit are stable in the charged state 0, -1 and 
-2. A comparison of the formation energies for oxygen 
and Zn vacancies reveal that the most stable defects are 
O vacancy in the Zn rich limit in the charge state -1-1 
(at p-type region) and also neutral O vacancies, while Zn 
vacancies in the O rich limit are found to be stable in 
the charge state -1 and -2 (at n-type region). It is in- 
teresting to note that for bulk ZnO, O vacancies in the 
charged state -1-1 exhibit negative U behavior and are 
never stable^. The reason for the negative U behav- 
ior is explained in Ref. [sO and we discuss below in some 
details. When one oxygen atom is removed from a per- 
fect bulk ZnO crystal, four Zn dangling bonds are created 
each contributing 1/2 electron to a neutral vacancy. This 
interaction results in a completely symmetric ai state ly- 
ing in the band gap and three almost degenerate higher 
energy states in the conduction band. For a neutral va- 
cancy, the ai state is occupied by two electrons and the 
energy is lowered when four Zn atoms surrounding the 
vacancy approach each other, resulting in 12% inward 
relaxation'^^. But, in case of -f 1 charge state of oxygen 
vacancy, this ai state is filled by one electron (half-filled), 
and the electronic energy gain is too small to overcome 
the strain energy, which leads to instability in the system. 
On the contrary, for the cluster considered in our work, 
the co-ordination of O is different from that of bulk. In 
the cluster one oxygen atom is coordinated with three Zn 
atoms leading to three dangling bonds upon creation of a 
vacancy. For a neutral vacancy, one defect state is occu- 
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pied by two electrons, which leads to large inward relax- 
ation (6.6%, 6.6%, 11.8%) similar to the bulk system^. 
For the +1 charged state, the defect state is half filled. In 
a non-spin polarized calculation the Fermi level is found 
to shift to a peak in a narrow structured density of states 
(due to the reduced co-ordination in the cluster) associ- 
ated with the defect. Such a high DOS at the Fermi 
level is conducive for the Stoner mechanism and favors 
spin-polarization to lower the overall bonding energy to 
make the system stable resulting in a magnetic moment 
of 1 fj,B- So the important difference between bulk ZnO 
and the cluster considered in this work lies in the fact 
that the defect state has a reasonable width (~ 0.9 eV) 
for bulk ZnO, which restricts the DOS at the Fermi level 
to a low value and does not allow gain in energy upon 
spin polarization. On the contrary, the defect states of 
the cluster with high value of DOS at the Fermi level 
offers gain in energy upon spin polarization^i. The low- 
ering of energy for oxygen vacancy in the charged state 
+1 is calculated to be ~ 130 meV upon spin-polarization 
concomitant with large relaxation (10. 7% outward, 10.7% 
outward, 6.9% inward) thereby adding to its stability. 
Similarly for the Zn vacancy the system becomes mag- 
netic in the charge state —1, and therefore becomes more 
stable. In contrast to the bulk, Zn vacancy in the cluster 
in the charged state —1 spans a larger energy range and 
the charged state —2 is only stable in the narrow n-type 
region, as shown in Fig. 9. So the magnetism is predom- 
inantly promoted by an unpaired electron in the charge 
state -|-1 and —1 for oxygen and Zn vacancies respectively 
and this mechanism may be responsible for magnetism in 
oxide nano-particles seen even in the absence of doping 
with magnetic elements^. 



V. CONCLUSION 



In conclusion, we have studied the energetics and mag- 
netic interactions in 3d transition metal doped ZnO 
cluster from ab initio calculations in the framework of 
GGA-I-U. Our calculations reveal that all the 3d transi- 
tion metal atoms couple antiferromagnetically in a pris- 
tine ZnO cluster with the possible exception of Ti. The 
presence of Zn and O vacancies have crucial effect on 
magnetic interaction. Some of the stable defect states 
(i.e. vacancy with some particular charge state) are 
found to stabilize ferromagnetism in the cluster when TM 
atoms are close to each other. Such defects also stabilize 
unusual ionic state of the dopant TM atom at the sur- 
face of the cluster in agreement with recent experimental 
results. A kinetic mechanism induced by spin conserved 
hopping is shown to mediate ferromagnetism in the sys- 
tem, where the unusual ionic state of dopant TM plays 
an important role. We also argue that defects in ZnO 
clusters even in the absence of TM doping may render 
the cluster magnetic. 
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